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ABSTRACT 

Deep Herschel PACS/SPIRE imaging and ^^C0(2-l) line luminosities from the IRAM Plateau de Bure Interferometer are combined 
for a sample of 17 galaxies at z > 1 from the GOODS-N field. The sample includes galaxies both on and above the main sequence 
(MS) traced by star-forming galaxies in the SFR-M* plane. The far-infrared data are used to derive dust masses, Mdust, following 
the Draine & Li (2007) models. Combined with an empirical prescription for the dependence of the gas-to-dust ratio on metallicity 
(^^gdrO^o))* the CO luminosities and Mdust values are used to derive for each galaxy the CO-to-H2 conversion factor, aco- Like in 
the local Universe, the value of aco is a factor of ~5 smaller in starbursts compared to normal star-forming galaxies (SFGs). We 
additionally uncover a relation between aco and dust temperature (rdust; Q^co decreasing with increasing Tdust) as obtained from 
modified blackbody fits to the far-infrared data. While the absolute normalization of the Qfco(7'dust) relation is uncertain, the global 
trend is robust against possible systematic biases in the determination of Mdust, ^^gdrO^o) or metallicity. Although we cannot formally 
distinguish between a step and a smooth evolution of ofco with the dust temperature, we can unambiguously conclude that in galaxies 
of near-solar metallicity, a critical value of Tdust = 30K can be used to determine whether the appropriate aco is closer to the "starburst" 
value (1.0 Mo(K km s'^ pc^)"!, when Tdust > 30 K) or closer to the Galactic value (4.35 M0(K km s"^ pc^)"!, when Tdust < 30 K). 
This indicator has the great advantage of being less subjective than visual morphological classifications of mergers/SFGs, which can 
be difficult at high z because of the clumpy nature of SFGs. Using Tdust to select the appropriate aco is also more indicative of ISM 
conditions than a fixed Ljr criterion. In the absence of far-infrared data, the off'set of a galaxy from the star formation main sequence 
(i.e., Alog(SSFR)Ms = log[SSFR(galaxy)/SSFRMs(M*,z)]) can be used to identify galaxies requiring the use of an ofco conversion 
factor lower than the Galactic value (i.e., when Alog(SSFR)Ms ^ 0.3 dex). 
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1. Introduction 

A correlation between the star formation rate (SFR) and the stel- 
lar mass (M*) of star-forming galaxies has been observed over 
the last lOGyr of lookback time (SFRcx , with 0.5 < a < 
1.0; "Brinchmann et al. 2004; Schiminovich et al.'2007VN oeske 
et al. 2007; Elbaz et al. 2007; Daddi et al. 2007; Pannella et al. 



2009 



2010 



Rodighiero et al.l|201Qt |Qliver et al.||2QlQt [Magdis et~ 



Karim et al.||2011[ [Mancini et al.||2011| ). The existence 



Send offprint requests to: B. Magnelli & A. Saintonge, e-mail: 
magnelli@mpe . mpg . de , amelie@mpe . mpg . de 

* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). Based also on observations carried out 
by the Herschel space observatory. Herschel is an ESA space obser- 
vatory with science instruments provided by European-led Principal 
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of this main sequence (MS) of star formation is usually inter- 
preted as evidence that the bulk of the galaxy population is form- 
ing stars gradually with long duty cycles. In this interpretation, 
galaxies situated on the MS are consistent with a secular mode 
of star formation, likely sustained by a c ontinuous gas accretion 
from the IGM and along the cosmic web ( Dekel et al. 2009} Dave 
|et al.||2Ql Q), while star-forming galaxies located far above the 
MS are strong starbursts with shortd uty-cycles, mainly triggered 
by major mergers (see, e.g. Engel et al. 2010 ). In the last years, 
independent observational tests have given weight to this inter- 
pretation. For example, it is observed that galaxies situated on 
the MS of star-formation have a disk-like morphology, relatively 
low SFR surface density (Wuyts e t al.„2Qllb ), cold dust tem- 
perature (Mag nelli et al.|2012[ see also, Magnelli et al. in prep.) 
and large polycyclic aromatic hydrocarbons (PAHs) equivalent 
width (EW; |Nordon et al.|2Q12| [Elbaz et al.|2QTT] ). On the other 
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hand, ofF-MS galaxies (i.e., situated above the MS) have cuspier 
morphology, higher SFR surface density, hotter dust temperature 
and smaller PAHs EW ( Wuyts et al. 201 lb; MagnelH et al. 2012 ; 
INordon et al.|2012t|Elbaz et al.|201 1| ). 

OfF-MS galaxies reaching very high specific star formation 
rates (SSFR=SFR/M*) do so by having enhanced star formation 
efficiencies (SFR/MH2; where is the molecular gas mass), 
as observed in local ULIRGs and high redshift submillimetre 
galaxies (e g. |Gao & Solomon|2004[ [Daddi et al.|2QT0[ [Genzel 



et al.||2010[ Saintonge et al.||2Qllb| ). These results are however 
complicated by the fact that the CO-to-H2 conversion factor. 
Of CO. also likely changes between normal star-forming di scs on 
the main sequence and m erger-driven starbursts ( Solomon et al. 
1997[ Tacconi et al.|2008[ ). Having either a measurement of aco 



in any given galaxy, or a reliable prescription to assign its value 
is therefore critical to study variations in star formation effi- 
ciency across the galaxy population. 

There are several methods to measure aco in high-redshift 
galaxies. For example a careful analysis of the stellar and dy- 
namical masses of luminous SMGs has been used to determine 
their CO-to-H2 conversion factor, aco (|Tacconi et al.||2008|). 
Although successful, this method requires very high quality in- 
terferometric CO observations (i.e., spectroscopically and spa- 
tially resolved) in order to determine the dynamical masses of 
galaxies. Consequently, before the advent in the near future 
of the full capabilities of the Atacama Large Millimeter Array 
(ALMA), this method will be only available for a handful of 
CO-bright galaxies. In order to constrain aco for a larger and 
fainter sample of galaxies another method is needed. In the local 
Universe, a method used to derive aco relies on the determi- 
nation of the gas mass of a galaxy via the measurement of its 
dust mass, using far-in frared observation s, and an assumed gas- 
to-dust ratio (see, e.g., Leroy et al.|2011| ). Thanks to the advent 
of the Herschel Space Observatory, this method can now be ap- 
plied toajarge sample of high-redshift galaxies. Magdis et aL] 
( 2011 20 1 2 ) have pioneered this method on a sample including 
a dozen high-redshift galaxies. Based on this recent progress we 
apply this method to a larger sample of high-redshift galaxies. 

Here we compile a large sample of 17 z > 1 galaxies in 
the GOODS-N field with both deep Herschel far-infrared mea- 
surements and accurate IRAM Plateau de Bure Interferometer 
(PdBI) line luminosities. After describing the sample selection, 
the data products, and the derived quantities in §Q4] we com- 
bine these measurements to measure the value of the conver- 
sion factor aco for each galaxy, discuss its trends with diff'er- 
ent global galaxy properties, and address the impact of possible 
systematic biases (Q. Finally, the origin and impact of the ob- 
served relations between aco and quantities such as Tdust and 
MS off'set are discussed in ^before the main conclusions are 
summarized in ^ 

All r e st-fram e quantities derived in this work assume a 



Chabrier 
r 



(2003) IMF, and a cosmology with //q 
Mpc-\Q^ = 0.27 and Da = 0.73. 



71km s 



2. Sample Selection 

Our goal is to study the link between molecular gas, dust, star- 
formation activity and stellar mass in high-redshift galaxies. 
Such a study requires using direct tracers of molecular gas and 
dust, and probing a large dynamic range in the SFR-M=, plane. 
One of the deep fields best suited for this study is GOODS-N, 
having a broad multi- wavelength coverage including some of the 
deepest Herschel maps (with both the R\CS and SPIRE instru- 
ment, see Section |3.1| ), and with its northern location, being a 



privileged target for the IRAM Plateau de Bure Interferometer 
(PdBI). The interferometric observations produce the CO line 
fluxes from which total molecular gas masses can be extrapo- 
lated, while the dust content of high-redshift galaxies can be 
studied using the wide far-infrared wavelength coverage from 
Herschel. 

We found 9 GOODS-N galaxies in the published litera- 
tur e with CO(l-O ) or C0(2-l) measurements, HDF169 (GN26) 
in IFray er et al.] ( [20Q8| ), BzK-4171, BzK- 21000, BzK-16000 , 
BzK-17999, BzK-12 591 and BzK - 25536 in [Daddi et al.1 ( [2QT0| ), 
HDF242 (GN19) in "Ivison et al.| ( [20TT] ), and finally GN20 in 
CariUi et al.]p010[ ). Since BzK-25536 is not detected in the deep 
PACS and SPIRE images (see Section [TT] ), it is excluded from 
our study. The remaining 8 galaxie^ are shown in the SFR-M* 
plane in the left panel of Figure[T](see Section [TT] for a complete 
description of the data used to derive SFR and M=,). All of these 
galaxies lie above the mean star formation main sequence, with 
three SMGs being the strongest outliers. In order to study the 
evolution of the molecular gas content of galaxies as a function 
of their location in the SFR-M* plane, we have obtained new 
IRAM PdBI CO data for galaxies located on or even below this 
main sequence of star-formation. 

At z~ 1, the Herschel observations of the GOODS-N field 
allow for the detection of galaxies with SFR> 10Moyr"\ 
probing well into the main se quence for galaxies with M,.> 
lO^Q- ^Me (see, e.g , Rodighiero et al.|2010| [Noeske et al.|2007 



Elb az et al.||2QQ7l |Karim et al.||2011|). Here we use the def - 



inition of the main-sequence from [Rodighiero et al.| pOlOj ), 



i.e., log(SSFR)Ms = a\og{M,)+l3 where (a , yS) = (-0.27, 2.6), 
(-0.51,5.3) and (-0.49,5.2) at 0.5 < z < 1.0, 1.0 < z < 1.5 
and z > 1.5, respectively. This choice is motivated by the use 
of the FIR as a star-formation indicator both in [Rodighiero et al.| 
( |2010 !) and in this study as well as by the use of fully consistent 
stellar mass estimateq^ However, none of our results strongly 
depend on this specific definition. 

We selected 9 IRAM PdBI CO targets from the com- 
plete volume-limited sample of galaxies with 1.0 < z < 
1.3, SFR> 15Moyr-^ and M,> 10iO-2Mo(see table [l]). In 
addition, we required the galaxies to have relia ble spectro- 
scopic redshift measurements from [Barger et al.[ ([2008 ). We 
note that our sample is not intended to be complete, but 
rather to sample as uniformly as possible the SFR-M=, plane 
(see Fig. [T]). The combined sample of 17 galaxies, includ- 
ing the objects extracted from the literature and our new 
CO targets, spans almost two orders of magnitude in SFR at 
fixed stellar mass; the final distribution of main sequence off'- 
set (i.e., Alog(SSFR)Ms = log[SSFR(galaxy)/SSFRMs(M„z)], 
where SSFR(galaxy) is the specific star- formation rate of our 
galaxies and SSFRms(M*,z) is the specific star-formation rate 
of MS galaxies at the redshifts and stellar masses of our galax- 
ies) is shown in the right panel of Fig.[T] 



^ We note amo ng these 8 galaxies, B zK-21000 and GN20 were ini- 
tially analyzed in Magdis et al. pOllj ), while the remaining 6 are also 



presented in Magdis et al. ( 2012F 

^ Results of Rodighiero et al. ( 2010 ) are co rrected for different IMF 
assumptions, namely from an Salpeter IMF in [Rodighiero et al.] ( |2010| ) 
to a Chabrier (2003 ) I MF in our study. Thus, SFRs and stellar masses 
from ^Rodighiero et al. ( [2010[ ) are divided by a factor 1.7 
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Galaxy 


Position 


Redshift 


Obs. CO 


^coAV 


Wco 


Flag^ 


^CO(l-O) 


Ref. CO 


log(M,) 


Metallicity'' 




RA 


DEC 




Transition 


Jy km s"^ 


km s"^ 




lO'^ Kkms-ipc2 




logCMo) 


12+log(0/H) 


PEPJ123712+621753 


12 37 12.15 


62 17 53.95 


1.249 


2-1 


0.15 + 0.04 


409 + 170 


1 


4.2+1.1 


This work 


10.62 


8.73 


PEPJ 123709+62 1507 


12 37 09.13 


62 15 07.92 


1.224 


2-1 


0.28 ± 0.06 


231 + 90 


2 


7.5 + 1.6 


This work 


10.36 


8.67 


PEPJ123759+621732 


12 37 59.47 


62 17 32.86 


1.084 


2-1 


0.44 + 0.11 


406 + 115 


1 


9.3 + 2.3 


This work 


10.58 


8.72 


PEPJ123721+621346 


12 37 21.45 


62 13 46.06 


1.021 


2-1 


0.39 + 0.13 


236 + 160 


2 


7.3 + 2.4 


This work 


10.53 


8.71 


PEPJ123615+621008 


12 36 15.19 


62 10 08.65 


1.027 


2-1 


< 0.42 




3 


7.9+1.5 


This work 


10.90 


8.78 


PEPJ123634+620627 


12 36 34.50 


62 06 27.94 


1.215 


2-1 


< 0.27 




3 


7.1 + 1.3 


This work 


11.06 


8.80 


PEPJ123633+621005 


12 36 33.65 


62 10 05.77 


1.016 


2-1 


0.78 + 0.19 


393 + 120 


2 


14.4 + 3.5 


This work 


10.97 


8.79 


PEPJ 123646+62 1141 


12 36 46.18 


62 11 41.99 


1.016 


2-1 


0.40 + 0.09 


157 + 60 




7.4+ 1.6 


This work 


11.06 


8.80 


PEPJ 123750+62 1600 


12 37 50.89 


62 16 00.69 


1.170 


2-1 


0.53 + 0.08 


649 + 1 10 




13.0+1.9 


This work 


11.31 


8.83 


HDF169 


12 36 34.53 


62 12 40.92 


1.224 


2-1 


3.45 + 0.93 


560 




92.3 + 24.9 


Frayer^ 2008] 


10.63 


8.73 


BzK-4171 


12 36 26.52 


62 08 35.30 


1.465 


2-1 


0.65 + 0.08 


530 




24.6 + 3.0 


Daddi4 2UiO 


10.72 


8.75 


BzK-21000 


12 37 10.64 


62 22 34.17 


1.523 


2-1 


0.64 + 0.07 


444 




26.0 + 2.8 


Daddi4 20i0 


10.99 


8.75 


BzK- 16000 


12 36 30.09 


62 14 27.96 


1.522 


2-1 


0.46 + 0.05 


194 




18.7 + 2.0 


Daddi4 2UiO 


10.98 


8.75 


BzK- 17999 


12 37 51.82 


62 15 20.11 


1.414 


2-1 


0.57 + 0.06 


440 




20.1+2.1 


Daddi4 20i0 


11.18 


8.82 


BzK-12591 


12 37 41.37 


62 12 51.18 


1.600 


2-1 


0.84 + 0.16 


600 




37.5 + 7.1 


Daddi4 20i0 


11.32 


8.79 


HDF242 


12 37 07.20 


62 14 08.08 


2.490 


1-0 


0.21+0.03 


850 




62.5 + 8.9 


Ivison^ZUiij 


11.41 


8.80 


GN20 


12 37 11.87 


62 22 12.46 


4.050 


2-1 


0.64 + 0.16 


670 




142.1 + 35.5 


Carini4 20ia' 


11.40 


8.80 



Notes. 

^ Significance of our PdBI detections: 1 - secure detection; 2 - tentative detection; 3 - upper limit. 

^ Metallicity are estimated using the mass-metallicity relation and converted into the Denicolo et al.|(2002) system. 
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Fig. 1. (Left) The z ~ 1 star formation main sequence from [Rodighiero et al . | ( [20 1 0| and the position of the GOODS-N sample 
with CO measurements and PACS/SPIRE photometry. Galaxies with new IRAM CO observations are shown as black circles while 
the red symbols are for galaxies with CO data from the literature d Frayer et al.|2QQ8 Daddi et al.|2010 Carilli et al.|201Q[ [Ivison 
et al.|2QlT] ). Among the red symbols, triangles stand for galaxies referred to SMGs. (Right) Distribution of distance from the star 
formation main sequence (AlogCSSFR)^^), for the galaxies with the new IRAM CO observations (dashed histogram), and the 
combined sample including galaxies CO measurements from the literature (open histogram). Bins of 0.25 dex in Alog(SSFR)M5 are 
used. 



3. Data 

3.1. Herschel and multi-wavelength observations 

The new IRAM PdBI CO targets were selected from the 
Herschel observations of the GOO DS-N field t aken as part of the 
PACS Evolutionary Probe (PEFp| |Lutz et al.| [2011) guaranteed 
time key programme. However, to characterize their FIR/submm 
SEDs we use here the Herschel observations of the GOODS- 
N field taken as part of the GOOD S-Herschel key programme 
(PI: D. Elbaz; E lbaz et al.||2Qll| ). This key programme pro- 
vides deeper PACS observations of the GOODS-N field along 
with deep SPIRE observations, allowing us to make use of the 
full wavelength coverage ( 100-500 yum) off'ered by the Herschel 



http://www.mpe.mpg.de/ir/Research/PEP 



Space Observatory. The GOODS -Herschel survey an d data re- 
duction methods are described in [Elbaz et al.| ( |20 lT] ). Here, we 
only summarise the properties relevant for our study. 

PACS and SPIRE source extraction was done using a point - 
spread- function (PSF) fitting analysis (Magnelli et al. 2009), 
guided using the position of sources detected in the deep 
GOODS-N 24 yum observations from the Multiband Imaging 
Photometer (MIPS; |Rieke et al.|2004|) onb oard the Spitzer Space 
Observatory (3cr ~ 20 /iJy; [Magnelli et al., 201 1) . This method 
both mitigates blending issues encountered in dense fields, and 
provides a straightforward association between MIPS, PACS 
and SPIRE sources as well as with the IRAC (Infrared Array 
Camera) sources from which the M IPS -24/70 yum catalogue was 
constructed ( MagneUi eFaL][20TT] ). The GOODS-N PACS and 
SPIRE catalogues reach a 3cr limit of 1.1, 2.7, 7.5, 9.0 and 12.0 
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Table 2. Mid- and far-infrared properties of our combined sample 



Galaxy 

PEPJ123712+621753 
PEPJ123709+621507 
PEPJ123759+621732 
PEPJ123721+621346 
PEPJ123615+621008 
PEPJ123634+620627 
PEPJ123633+621005 
PEPJ123646+621141 
PEPJ 123750+62 1600 
HDF169 
BzK-4171 
BzK-21000 
BzK- 16000 
BzK- 17999 
BzK-12591 
HDF242 
GN20 



^24 

153 ±5 
128 ±6 
217 ±7 
238 ±8 
327 ±7 
108 ±9 
572 ±7 
300 ±6 
188 ±6 
445 ±7 
140 ±6 
382 ±6 
190 ±6 
229 ±8 
377 ±6 
255 ±9 
68 ±4 



^70 

mJy 

<2:o 

< 2.0 

< 2.4 

1.8 ±0.6 
4.5 ± 0.6 

< 2.0 
2.7 ± 0.6 

< 2.0 

< 2.0 
13.2 ±0.7 

< 3.2 

4.9 ± 0.7 

< 2.0 

< 2.0 

< 3.3 

< 2.0 

< 2.0 



100 

mJy 

2.6 ±0.3 

2.1 ±0.4 

1.2 ±0.4 
3.1 ±0.3 

4.3 ± 0.4 
7.1 ±0.4 

8.5 ±0.3 
2.9 ±0.3 

1.7 ±0.3 
34.2 ± 0.6 
2.9 ±0.3 
8.1 ±0.6 

1.6 ±0.6 
4.1 ±0.5 

9.8 ± 0.4 
1.1 ±0.3 
0.7 ± 0.4 



^160 

mJy 
3.5 ±0.7 
3.0 ± 1.3 

5.3 ± 1.7 
8.5 ±0.8 
9.2 ± 1.0 
13.4 ± 1.2 
16.3 ± 1.0 
5.8 ±0.7 
4.8 ± 1.1 

63.2 ± 1.6 
9.0 ± 1.0 
15.1 ± 1.4 
3.7 ±0.9 
10.8 ± 1.4 

18.3 ± 1.1 

6.4 ± 1.3 
5.4 ± 1.0 



^250 

mJy 
6.2 ± 2.5 
8.0 ±2.5 

10.0 ± 2.5 
14.5 ±2.5 
15.9 ±3.1 

18.1 ±4.0 

20.2 ± 2.5 
12.7 ±2.5 
8.6 ±2.5 
62.5 ± 3.7 

14.5 ± 2.5 
24.4 ± 2.5 
9.2 ± 2.5 
14.7 ± 2.5 
24.9 ± 2.5 
23.9 ±4.6 

18.6 ±2.7 



^350 

mJy 



< 

10.8 
14.3 
16.1 
13.0 
10.7 

8.8: 

47.1 

< 

20.1 

< 

11.2 
20.5 
28.3 
41.3 



9.0 
9.0 
9.0 
±3.0 
±3.0 
±3.0 
±3.0 
±3.0 
±3.0 
±3.0 
9.0 
±4.7 
9.0 
±3.0 
±3.0 
±6.4 
±5.2 



"^500 

mJy 

<T2:o 

< 12.0 

< 12.0 
8.0 ±5.0 

< 12.0 
15.1 ±4.0 

< 12.0 

< 12.0 

< 12.0 
19.6 ±4.0 

< 12.0 
11.6 ±7.4 

< 12.0 

< 12.0 
9.9 ± 4.0 
23.8 ±4.0 



"^850 

mJy 



2.2 ±0.8 



8.0 ±3.1 



1160 

mJy 

< 1.7 

< 1.6 

< 2.1 

< 1.6 

< 1.7 

< 2.0 

< 1.6 

< 1.7 

< 1.8 

< 1.8 

< 1.7 

< 2.2 

< 1.7 

< 1.8 

< 1.7 

< 1.6 



"^2200 

mJy 



"^ 3300 

mJy 



0.32 ±0.15 0.04 ±0.06 



39.7 ±6.1 20.3 ±2.0 10.5 ± 0.7 0.9 ±0.15 0.33 ± 0.06 0.027 ± 0.( 



Notes. 

Sxx correspond to flux densities at the observed wavelength XX in microns. 
^ Upper limits corresponds to 3cr upper limits. 
^ These flux densities are taken from Pope et al. f2006|. 
^ These flux densities are taken from Penner et al. ( 2 01 It. 

These flux densities are taken from |Magdis et al.,^2011) . 



mJy at 100, 160, 250, 350 and 500 yum, respectively. 

The GOODS -N field benefits from an extensive multi- 
wavelength coverage. A PSF-matched catalogue was created for 
GOODS-N as part of the PEP surve}0( Berta et al. 2010, 2011), 
including photometry from far-UV to near-infrar ed and a coUec - 
tion of spectroscopic redshifts (mainly from Bar ger et al.|2 008^, 
complemented when needed with photometric redshifts com- 
puted using EAZY (Brammer et al. 2008 ). Stellar masses were 
derived by fitting Bruzual & Chariot ( 2003 ) templates to these 
far-UV to near-infrared photometries (see Section |4.3| for more 
detail). We cross-matched our IRAC-MIPS-PACS-SPIRE cata- 
logue to this multi-wavelength catalogue, using the IRAC posi- 
tions of each catalogue and a matching radius of T'. 

Galaxies with CO observations from the literature were also 
cross-matched with the UV- to-far-infrared multi- wavelength 
catalogue. We verified that the available photometry and spec- 
troscopic redshift from the literature were consistent with those 
from our multi- wavelength catalogue. We find no systematic dis- 
crepancies. Finally, we cross-matched using a 9''matching radius 
our sample with the 1.16 mm catalogue of [Penner et al.|p011 



which was built by combining the 1 . 1 mm AzTEC ([Perera et a 



2008) ) and the 1.2 mm MAMBO (Greve et al.||2008| maps. We 



found a clear 1.16 mm detection for GN20 and obtained mean- 
ingful 1.16 mm upper limits (i.e., 3cr flux limits at the location of 
our galaxies) for the rest of our sample. Table|2]gives the infrared 
photometry for all the galaxies in the combined sample. 

3.2. IRAM PdBI observations 

Observations of the 9 //^r^c/z^/- selected nor mal star-forming 
galaxies were performed at the IRAM PdBI ( [Guilloteau et al.| 
[1992^ in 2010 July-October under average summer conditions. 
We observed the C0(2-l) line (rest frequency of 230.538 GHz), 
which is redshifted into the 3mm observing band for sources at 
z ~ 1 . At the time of observing, the PdBI array was in compact 
configuration with 5 antennae available. The WideX correlator 
was used, providing dual polarization data with spectral resolu- 
tion of 1.95 MHz over a total bandwidth of 3.6 GHz. 



The data were calibrated with CLIC within the standard 
IRAM software package, GILDAS|^ A standard passband cal- 
ibration scheme was applied, followed by phase and amplitude 
calibration. The flux calibration was done using observation of 
reference bright quasars. Due to the poor observing conditions 
during some of the runs (e.g. high precipitable water vapour, 
strong winds, low elevation), particular care was taken to flag 
data with high phase noise. Finally, data cubes were created us- 
ing the MAPPING environment and cleaned using the CLARK 
version of CLEAN implemented in GILDAS. The typical RMS 
in the data cubes is of 0.8 mJy/beam with a typical beam di- 
ameter of -5'\ The data cubes were examined for sources at 
the expected spatial and spectral positions. The integrated line 
map and spectrum of each target are shown in Figures [2] and 
[3] Four of the sources are securely detected, three are tentative 
detections, and the remaining two are non-detections. For the 
detections, line fluxes were measured by integrating the spectra 
over the width of the signal (Wco)- Due to the absence of con- 
tinuum in our spectra, Wqo is defined as the frequency range 
maximizing the signal to noise ratio (S/N) (i.e., the peak of the 
S/N versus Wqo plot). This Wco also corresponds to the fre- 
quency range where the curve of growth of the line flux reaches 
its plateau (see center right and right panels of Figure [2] and [3]). 
To create the S/N versus Wqo and line flux versus Wqo diagrams, 
we have to fix the central frequency around which we increase 
Wco- These fixed frequencies are defined as being the central 
frequencies of the yellow areas displayed in the left panels of 
Fig.[2]and|3] Uncertainties on Wqo are defined as being the range 
of Wco over which S/N is greater than max[S/N(Tyco)]-l- For 
non-detections we measured 3cr upper limits as suming Wco = 
300kms"\ typical for high-redshift galaxies ( [Tacconi et al. 
[2010 ). A summary of the observed CO fluxes and Wco for both 
literature and new sources is presented in Table [T] 



publicly available at http://www.mpe.mpg.de/ir/Research/PEP/ 



http://www.iram.fr/IRAMFR/GILDAS 
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4. Analysis 

4.1. The CO(1-0) luminosities 

Since the standard techniques to extrapolate from CO mea- 
surements are calibrated for observations of the / = 1 ^0 
transition, the fluxes coming from observations done in higher 
transitions need to be converted into fluxes of the / = 1 ^0 
transition. All but one of the galaxies have observations of the 
J = 2 ^ I transition, with HDF242 (aka GN19) observed di- 
rectly in / = 1 ^ transition pvison et al.||2011| ). Thanks to 
this homogeneity, our combined sample is mostly free of un- 
certainties in the excitation corrections when combining multi- 
transition samples. We adopt a ratio of S (2-i) / S (i-o) = 3, which 
is typically observed within the disks o f local normal star- 
forming galaxies (e.g., L eroy et al.||2009 ) as well as in normal 
high-redshift st ar-forming galaxi es (i.e., BzK gala xies; see e.g., 



Aravena et al 



2Q10f l Recently 



BothweU et al. 



( 2012] ) found 



5(2-i)/5'(i_o) = 3.36 ± 0.52 for a sample of 40 luminous SMGs, 
i.e. galaxies mostly situated far above the MS of star-formation. 
Because there is a clear agreement between the S (2-1) /S (i_o) flux 
ratio found in normal SFGs and SMGs, we conclude that the 
conversion of CO / = 2 ^ 1 fluxes into CO / = 1 ^ fluxes 
does not introduce any difl'erential biases in the CO(l-O) lumi- 
nosity estimates of on- and ofl'-MS galaxies. The CO(l-O) lumi- 
nosities can then be calculated following [Solomon et al. ( 1997 ): 

Vr^r.n_n. = 3.25 X 10^^(i_o)Ay Vn%,A'(l + z)"', (1) 



'CO(l-O) 



where 5'(i_o)Ay is the integrated line flux in Jykms"\ 
y(i_o) = 115.271 GHz is the rest- frame frequency of the 
CO(l-O) line, and Z)l is the luminosity distance in Mpc. Table [T] 
summarises the inferred CO properties of our sample. 



4.2. Dust masses, dust temperatures and SFRs 

To infer the dust ma ss of our galaxies we use the dust model of 
[Draine & Li| ( [2007] hereafter DL07). This model describes the 
interstellar dust as a mixture of carbonaceous and amorphous 
silicate grains with a size distribution mimicking the Milky Way 
(MW) extinction curve. The smallest carbonaceous grains have 
the properties of PAHs and their abundance is parametrized 
by the RAH index, ^pah, i.e., the fraction of dust mass in the 
form of PAH grains with less than 10^ carbon atoms. The DL07 
model assumes that a large fraction of the dust (i.e., I - y with 
0.0 < y < 0.3) is situated in the difl'use ISM and exposed to a 
radiation field with a constant intensity f/min (U being a dimen- 
sionless scale factor normalized to the local interstellar radiation 
field of the MW). The remaining fraction of the dust (i.e., y) is 
exposed to a power-law (a) radiation field ranging from t/min 
to t/max and corresponds to the dust enclosed in photodissocia- 
tion regions (PDRs). Combining these two components, the dust 
emission spectrum of a galaxy is given by : 



pectrum or a galaxy is given oy : 

[(^-7)Py\(lFAH, t/min) + yPviqFAU, U^in, ^/max, Q^)](2) 

where p^y\qmn, ^min) is the emitted spectrum per unit dust mass 
for the component heated by a single radiation field t/min, and 
Pviqmn, Ujrnn, ^^max, Of) is the cmission spectrum for the PDR 
dust. DL07 found that the exact value of a does not have a large 
influence on the quality of their fits and fix a = 2. In addition. 



^ It corresponds to a ^co(2-i)/^co(i-o) = ^^(2 l)/Tb(l 0) ratio 
of 0.75, where Tb is the equivalent Rayleigh- Jeans brightness tempera- 
ture. 



they suggest to fix Umax = 10^ and to use 0.7 < Uj^nn ^ 25 when 
submillimeter data are not available. This lower cutoff' for Uj^in 
is made to avoid erroneously fitting very cold dust components 
to a set of data with no submm constraints. DL07 conclude that 
this cutoff' does not introduce any systematic biases in their dust 
mass estimates. 

Following these prescriptions, we build a grid of models with 
diff'erent PAH abundances (0.47% < ^pah < 4.6%), values of 
^min (0.7 - 25) and values of y (0.0 - 0.3). At each grid point, we 
compare the model SED and the FIR/submm photometry, with 
the dust mass given by the normalization of the SED that min- 
imizes;^^. Integrating this best-fitting SED template from rest- 
frame 8 to lOOOyum, we also infer for the grid point the infrared 
luminosity (Ljr) and SFR (with SFR[Moyr-i] = lO'^^ Lir[L0], 
assuming a Chabrier IMF). We assign for each galaxy a value 
of Mdust. ^iR and SFR given by the geometric mean of the re- 
spective quantity across the region of the parameter space where 

< ~^ ^' These quantities, along with their Icr uncertainties 
derived over the same grid points, are given in Table [3] 

We find that our grid of DL07 models provides a good fit 
to the photometric data points, with a median x^ = 1.8 for 
A/dof = 1. This is in agreeme nt with DL07, [Munoz-Mateos et al 
( 2009') and ' Dale et aLl ( |2012 ), who found that this particular set 
of models could reproduce well the large wealth of data avail- 
able for the SINGS and KINGFISH galax ies. We note t hat some 
of our galaxies are also included in (Magdis et al. 2012 ; i.e., the 
BzK galaxies). Our dust mass estimates for these galaxies are in 
agreement with those of |Magdis et al.| ( [2QT2| ) who also used the 
DL07 model. 

Among our 17 galaxies, 12 have clear (S/N> 3, 8 galaxies) or 
tentative detections (3 > S/N > 2, 4 galaxies) of Herschel ffuxes 
up to the rest- frame 160yum wavelength. For these galaxies our 
dust mass estimates should be robust. Indeed, iDale et al.| ( |2012| ) 
found that the lack of /Irest > 160 yum observations does not intro- 
duce any systematic biases, although the scatter in the resulting 
dust masses may increase by up to ~ 50%. The main conclusions 
of this paper are based on the comparison of our dust and gas 
measurements with empirical relations presented in'Leroy et al.| 
( 201 1 ), where results are also based on FIR measurements short- 
ward of rest- frame 160/im. This consistency with the method- 
ology of Ler oy et al.| ( |2011 ) trumps possible worries regarding 
the overall normalization of our dust masses in the absence of 
submm data. 

For the remaining 5 galaxies, reliable Herschel ffuxes are 
only available up to rest-frame lOOyum. To test for a systematic 
bias in the dust mass estimates of these galaxies, we re-fftted the 
other 12 galaxies this time omitting any ffux measurements at 
/Irest > lOOyum. Wc ffnd no systematic biases, although the scat- 
ter in the resulting dust masses increases by ~ 50% compared 
to the reference values obtained using the ffux measurements at 

/Irest > lOOyUm. 

In order to obtain a proxy of the luminosity-weighted dust 
temperature of each galaxy, we fft their far-infrared SEDs, in the 
optically thin approximation, with a single modiffed blackbody 
function: 



Sy OC 



Qxp(hy/kT^^,t) - 1 ' 



(3) 



where 5'y is the ffux density, 13 is the dust emissivity spectral in- 
dex and Tdust is the dust temperature. We ffx the dust emissivity 
spectral index to 1.5, a standard value shown to provide a good 
fit to the far-infrared SEDs of high-redshift galaxies (see, e.g., 
Magnelli et al.|2012| ). Since a single temperature model cannot 
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Fig. 2. C0(2-l) spectrum (left), integrated line map (center left panel), line flux versus Wco (center right panel) and signal over noise 
(S/N) versus Wco (right ) of our 7 PdBI targets with secure (PEPJ123712+621753, PEPJ123759+621732, PEPJ123646+621141 
and PEPJ123750+621600) and tentative line detection (PEPJ123709+621507, PEPJ123721+621346 and PEPJ123633+621005). In 
each spectrum, the yellow area displays the part of the spectrum over which the line fluxes are integrated. Dashed lines in the right 
two panels outline this value of Wco- To create the center right and right panels, we fixed the central frequency when increasing Wco 
to the central frequency of the yellow area displayed in the left panel. Dotted-dashed lines and triple-dotted-dashed lines represent 
the line flux and error of each object, respectively. Beam sizes are displayed in the lower left corner of each integrated line map. 
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Table 3. Dust properties of our combined sample 



Galaxy 


log(MjL07) 






j-BB 
dust 




Mo 


Lo 


Mo 


K 



PEPJ123712+621753 


8.91 


+ 


0.42 


11.64 


±0.05 


7.93 


±0.57 


31.5 


± 


6.5 


PEPJ123709+621507 


8.96 


± 


0.40 


11.58 


±0.06 


8.50 


±0.60 


25.0 


± 


5.0 


PEPJ123759+621732 


8.94 


± 


0.27 


11.46 


±0.05 


8.27 


±0.28 


27.0 


± 


3.0 


PEPJ123721+621346 


9.07 


+ 


0.11 


11.65 


±0.01 


8.45 


±0.10 


26.0 


± 


1.0 


PEPJ123615+621008 


9.33 


+ 


0.12 


11.78 


±0.01 


8.50 


±0.10 


26.0 


± 


1.0 


PEPJ123634+620627 


8.56 


+ 


0.16 


11.95 


±0.03 


8.55 


±0.10 


29.0 


± 


1.0 


PEPJ123633+621005 


8.61 


± 


0.05 


11.87 


±0.01 


8.32 


±0.07 


30.0 


± 


1.0 


PEPJ123646+621141 


9.26 


± 


0.15 


11.61 


±0.01 


8.52 


±0.12 


24.0 


± 


1.0 


PEPJ 123750+62 1600 


8.97 


+ 


0.25 


11.58 


±0.03 


8.35 


±0.25 


26.5 


± 


2.5 


HDF169 


9.03 


+ 


0.05 


12.61 


±0.01 


8.60 


±0.05 


37.0 


± 


0.5 


BzK-4171 


8.65 


± 


0.12 


11.97 


±0.03 


8.18 


±0.18 


34.5 


± 


2.5 


BzK-21000 


9.30 


± 


0.02 


12.38 


±0.01 


8.52 


±0.03 


33.5 


± 


0.5 


BzK- 16000 


8.96 


± 


0.32 


11.84 


±0.04 


8.43 


±0.43 


29.0 


± 


4.0 


BzK-17999 


8.67 


+ 


0.11 


12.03 


±0.02 


8.25 


±0.15 


34.0 


± 


2.0 


BzK-12591 


8.83 


+ 


0.05 


12.42 


±0.01 


8.42 


±0.03 


36.5 


± 


0.5 


HDF242 


9.30 


± 


0.06 


12.65 


±0.01 


9.02 


±0.12 


33.0 


± 


2.0 


GN20 


9.66 


± 


0.01 


13.23 


±0.01 


9.15 


±0.05 


41.0 


± 


0.5 



fully describe the Wien side of the far-infrared SED of galax- 
ies, because short wavelength observations are dominated by a 
warmer or transiently heated dust component, we exclude from 
our fitting procedure flux measurements with A^Qst < 50 jim. 

The absolute dust temperature derived using a single black- 
body model strongly depends on the assumed dust emissivity. 
For example, using JS = 1.5 leads to dust temperatures system- 
atically higher than if using j3 = 2.0 (AJdust ~ 4 K). Nonetheless, 
fixing 13 to 1.5 or allowing it to vary freely between 1.5 - 2.0 
does not significantly aff'ect the main results presented in Section 
|5] In any case, to prevent our results from only relying on the 
dust temperature, we also use the far-infrared color of galaxies 
(i.e., /y[60yum]//y[100yum]; hereafter /60//100) as an indepen- 
dent measure of the conditions prevailing in their ISMs (see, e.g.. 
Dale et al.|2001| ). Here, /60//100 is simply obtained using a log- 



linear interpolation between all our observed far-infrared data 
points. Because our observations always encompass the rest- 
frame 60 and 100 yum passbands, these estimates are robust and 
model independent. 



From the normalization of the modified blackbody function 
we can also compute a dust mass: 



(1 + z)K(yo)By(Ar:QsU ^dust) 



vo 



(4) 



where K(yo) is the dust mass absorption cross section at the ref- 
erence frequency vq. To allow direct comparison with the dust 
masses estimated using the DL07 model, we use the same dust 
mass absorpti on cross sect ion at vq = c/lOOyum, i.e., /c(yo) = 
3.13 m^kg"^ ( Li & Draine,,200 1 ). Dust temperatures and dust 
masses inferred from this single modified blackbody function 
are given in Table [3] Here again we quote the geometric mean 
across the region of parameter space satisfying < x^i^ + 1 and 
from which we define the Icr uncertainties. 

Fig. |4] compares the dust masses obtained using a single 
modified blackbody function to those inferred using the DL07 
model. Estimates based on a single modified blackbody function 
are systematically lower by a factor ~ 3. This finding is in perfect 
agreement with that of Dal e et al.| (|2012 ) based on local galax- 
ies observed as part of t he SINGS and KINGFISH surveys (see 
also |Magdis et al.|2012| ). These discrepancies arise from the fact 
that single temperature models do not take into account warmer 
dust emitting at shorter wavelength. The behavior is also con- 
sistent with the analytic radiative transfer work of [Chakrabarti] 
,& McKee ( 2008| ), who found that in high-redshift star-forming 
galaxies, dust masses inferred from isothermal models (equiv- 
alently single temperature blackbody approximation) are lower 
by a factor of ~ 2 than those inferred using multi-r components 
(see equation A8 of Chakrabarti and McKee 2008 ). In the rest 
of the paper we work with dust masses inferred using the DL07 
model as they account for the warmer dust emitting at shorter 
wavelength. This choice is also motivated by the need to reduce 
any possible correlations between our dust temperature and dust 
mass estimates (see Section|5^for more details). 
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Fig. 4. Comparison of the dust masses estimated using the DL07 
model with those estimated using a single temperature modified 
blackbody function. The solid line corresponds to the one-to-one 
relation while the dashed line is off'set by a factor 3. We quote 
the geometric mean of Mdust across the region of parameter space 
whQTQx^ < xliin + 1 f^^^ which we define its Icr uncertainty. 



4.3. Stellar masses and metallicities 

In order to estimate the stellar masses of our galaxies we use 
the multi- wavelength catalogue built by the PEP consortium and 
presented in Berta et al. ( 2010, 2011 ). It is a z + ^-selected cat- 
alogue containing photometry in 16 bands from the Galex far- 
UV to Spitzer near-IR wavelengths. Stellar masses were calcu- 
lated by fitting these multi-wavelength photometry to Bruzual &] 
[Charlot| ([2003 ) templates using FAST (Fitting and Assessment 
of Synthe tic Templates; Kriek et al. 2009^ an d prescriptions pre- 
sented by [Wuyts et'aL] ( |201 lb| ). We limit the [Bruzual & Chariot 



2003| ) templates to models with exponentially declining SFHs 



and a minimum ^-folding time of 300 Myr. This allows for a 
much better agreement between SFRs derived from optical-to- 
near-IR SED fits and those derived using mid/far-infrared obser- 
vations. Full details are presented in |Wuyts et al.| ( |201 la| ) and 



Wuyts etal. (2011b). 



Metallicities are derived using the stellar-mas s- metall icity 
relation at the redshift of the source (Genzel et al. [2012} see 
also E rb et al.|2006t P. Buschkamp et al., in prep.; Shapley et al. 
[20051 Liu et aL ]2008j : 



Zo = 2.181og(M,) - 0.0896 log(M,r + a. 



(5) 



where a = -4.45 if 1.0 < z < 1.5 and a = -4.51 if z > 1.5 
All metallicities are converted into the |Denicol6 et al.| ([2002|) 
calibration system using the conversion function of Kewley &[ 
[Ellison (2008 ). The Denicolo et al. calibration system provides 
the best agreement, at high stellar masses, between all different 
metallicity calibrators ( [Kewley & Ellison[2008[ ). Table[3]presents 
our inferred metallicities. We note that metallicities estimated 



using the "fundamental metallicity relatio n" of Mannucci et al. 
( |2010| ) extended to z ~ 2 in [Genzel et al.| ( |2012| ) are consistent 
with those derived using equation [5| 



5. The CO H2 conversion factor 

The total molecular gas mass is commonly derived from the 
measured CO line luminosity using a conversion factor ofco- 



= o:coLr 



'CO(l-O)' 



(6) 



where L' 



'CO(l-O) 



is in (K km s ^ pc^). The exact value of ^co 
and its possible dependency on a range of parameters were and 
are currently being actively investigated, both with observations 
and simulations (e.g., [Solomon e t al. 1987; Dickman et al. 19861 
Maloney & Black; 1 988 [|IsraeIl2000,,Bolatto et al.,2008HLeroy 
et al. 11201 1 i Mag dis et al.[|201H [Genzel et al.[|2012| [Glover & 



Mac Low 20 11nN^ayanan et al.[[2011t [Shetty et al.[[2011a|bT 



Schruba et al.i2012| ). 



There is a general consensus that within the MW ar n= 4.34 



Mn(K km s PC ) , with little spatial varia tions (e^. [Strong 



& Mattox|l996l[Dame et al.|2001 



Abdo etal. 2010r There is 



also ample evidence for at least two sets of circumstances where 
a CO departs significantly from the Galactic value: low metallic- 
ity environments and the very dense central regions of merg- 
ing systems. In low metallicity conditions, the CO molecules are 
lacking the dust shielding they require to survive the ambient ra- 
diation field everywhere but in the cores of GMCs, while the H2 
molecules can survive in self-shielded envelopes around these 
cores (e.g. [Wolfire et al.[[201()| ). The result of this preferential 
photodissociation of the CO molecule are values of aco larger 
by a factor of at least 5-10 (or even more) in local low-metallicity 
galaxies such as the SMC (Israel 1997; Joha nsson et al.[[1998[ 
Leroy et al.|2QTT|[Bolattoet al. 2011; Schrub a et al.|2012[ ). On 
the other hand, in systems undergoing strong starbursts triggered 
by major mergers, such as nearby ULIRGs and high-z SMGs, 
the value of aco is reduce d by a factor of 5 (Solomon et al. 
1997; Downes & Solomon [1958) Tacconi et al. [20Q8t but see 
Papadopoulos et al. 2012). 

While these boundary conditions are generally accepted and 
routinely applied, the transitions between these diff'erent regimes 
are still being explored. Here we combine the IRAM molecular 
gas data and the Herschel dust masses to estimate a conversion 
factor for each galaxy in our combined sample, and compare the 
values with different prescriptions for aco variations. 

5.1. The gas-to-dust ratio method 

Assuming that we can predict the gas-to-dust mass ratio ^gdr 
in a galaxy, and having accurate molecular gas and dust masses, 
the conversion factor a^co can be estimated as: 



o:co 



^GDR(/^o)-^dust 
^CO(l-O) 



(7) 



To es timate ^gdr for each galaxy, we use the relation of [Leroy 
et al., ( 2011) which relates it to metallicity, i.e., yUo: 

log ( WOt/o) ) = -0.85 X yuo + 9.4. (8) 



Recently, [Magdis et al.[ ( [201 1| ) have shown that this "dust- to- 
gas ratio method" constrains aco in two well-studied GOODS- 
N galaxies to values consistent with standard expectations. 
However, we wish to point out three important assumptions re- 
quired to extend this method from the local Universe to our high 
redshift galaxies: 



^ We note that all the values of aco quoted in this paper already in- 
clude an upward factor of 1.36 to account for the presence of Helium 
coexisting with the molecular hydrogen. 
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Since the CO observations were done in large part with the 
PdBI in compact configuration, and since the Herschel PSF 
even at the shortest wavelengths has FWHM- 6''(i.q., corre- 
sponding to linear scales of ~ 48 kpc at z ~ 1), we can only 
work with integrated dust and gas measurements. Therefore, 
we must rely on the assumption that the total C0(2-l) line 
fluxes and the FIR fluxes we measure are emitted from the 
same physical regions. In the local Universe, IR emission 
and CO measurements are seen to correlate well, with a tight 
relation between Lfir and Leo (e .g. 'Sanders & Mi rabel 1985[ 
[Sanders et aL][T99Tl |Gao & S olomon 2004 ), which can be 
interpreted as evidence for them being produced in the same 
region. Furthermore, in local normal star-forming galaxies, 
the radial profiles of C0(2-l) line emission and of IR flux. 



2. 



5.2. Results 

Figure [5] shows the value of ^co as a function of metallicity 
(yuo), dust temperature (rdust), far-infrared color (/60//100) and 
distance of a galaxy with respect to the MS of star-formation 
(Alog(SSFR)M5). Our galaxies roughly follow the aco-J^o rela- 
tion found in [Genzel et al. | ( [20T2 ^) , although with a larger scatter 
(~ 0.5 dex as opposed to the expected ~ 0.2 dex) and over a lim- 
ited metallicity range which makes the slope uncertain. We note 
that this larger scatter is expected because our sample contains 



normal SFGs and merger-driven starbursts as opposed to |Genzel 
|et al.| ( |2012| ) who studied normal SFGs. This scatter however 
correlates with r^ust^ /60//100 and Alog(SSFR)M5 . We note that 



as traced b y 24yum emission, are strikingly similar ( [Leroy 
let al.|2008|) . 

Equation[7|is built under the assumption that » Mm, and 
therefore that Mgas - Mii^ = ^gdrO^o) ^dust- While in nearby 
galaxies, on average Mhi~ SMh^, albeit wit h large galaxy- 
to-galaxy variations (' Saintonge et al.|2011a| ), there are mul- 
tiple reasons to suspect that the balance between atomic 
and molecular gas is significantly diff'erent in high redshift 
galaxies (see e.g., Obreschkow & Rawlings 2009| ). First, the 
molecular gas mass fraction = Mii^/(M^ -\- M^^), ob- 
served locally to be on average 6% in star-forming galax- 
ies ( [Saintonge et al. 201 1a), is found to be 30 - 50 % at 
z = 1 - 2 (e .g., |Tacconi efa l. 2010; Daddi et a l. 2010; Geach| 
et al.pOlT] ). It is therefore difficult to accommodate a sig- 



nificant atomic gas component within the total mass budget 
of these high-z galaxies as inferred from dynamical analy- 
sis. Observations of damped Lyman-a systems also suggest 
that the evolution of Qhi b etween z = and z = 2 is not a s 
strong (see e.g., figure 2 of Obreschkow & Rawlings|20()9) . 
Finally, both models and observations indicate that above a 
critical surface density of ~ 10 M© pc"^ most of the gas is 
found in molecular form, with HI saturating at that threshold 
(e.g.,|Wong & Blitz|2002[|Schaye|2004l|Blitz & Rosolowsky] 



2006[|Bigiel et al.|200 8V Since high-z star-forming galaxies 
have high surface densities, we can therefore expect minimal 
contributions from HI over the regions of the disks where the 
CO and FIR emission are originating. 
We also work under the assumption that the local relation 
between ^gdr and metallicity applies in the denser ISM of 
high-redshift galaxies. In denser environments, the deple- 
tion of h eavy element s from the gas phase is expected to 
increase ( [Jenkins |[2009^ ), leading to lower ^gdrO^o) values. 
However, in the denser environment of local-ULIRGs, only 
a mild decrease of Sgdr(j^ o) is observed, compared to nor- 
mal star- forming galaxies (jSolomon et al.|1997{ Wilson etaL] 
[2008; Santini et al. 20f0| [Clements et al.|2010| ). Thus, even 
if high-redshift galaxies exhibit high surface densities, the 
eff'ect on ^gdr(/^o) should not exceed a factor 2. The conse- 
quence of this eff'ect is further discussed in Section [5^ Due 



to their higher surface densities, high redshift galaxies might 
also have higher dust emissivity (see, e.g.,[Arce & Goodman 



T999l [Dutra et al.|2003||Cambresy et al.|2005||Michatowski 



et al. 2010| ) than local galaxies. This would affect our dust 
mass estimates. However, we show in Section \5?2\ than this 
effect does not strongly aff'ect our main results. 



Based on these assumptions, we can combine the IRAM molec- 
ular gas data and the Herschel dust masses to estimate aco in 
each of the GOODS-N galaxies. 



Gracia-Carpio et al7('2011') also proposed, in the local Universe, 
for a smooth evolution of a^co with feo/fioo- Their empirical 
relation is shown in the o^co-/6o//ioo panel of Fig. [5] Despite 
some possible discrepancies in the absolute normalization of our 
two correlations, we find that the relative variation of aco with 
/60 //loo found in our two studies are in excellent agreement. Our 
aco - Alog(S SFR)m^ relation is also consistent with the results of 
[Magdis et"aL]p012[ ) who find that galaxies situated off' the MS 
exhibit lower aco than galaxies on the MS of star-formation. 

In order to characterize the correlation of aco with the dif- 
ferent quantities, we fit these relations with (i) a constant, (ii) a 
step and (iii) a linear function. To take into account errors both 
in X and Y, we use a Monte-Carlo approach. We create 1 000 
mock samples with the same number of sources as our original 
sample. For each realization and for each galaxy, we attribute 
new values of aco, ^dust, /60//100 and Alog(SSFR)M5 selected 
into a Gaussian distribution centered at their original values and 
with a dispersion given by their errors. We then fit each Monte- 
Carlo realization with our three diff'erent models (i.e., constant, 
step and linear functions). Table [4] summarizes the results. For 
each model, we quote the mean value and dispersion of each 
fitting parameter across the 1 000 Monte-Carlo realizations. To 
assess the degree of correlation between aco and Tdust, /60//100 
and Alog(SSFR)M5 , we also perform Spearman's test on each of 
the Monte-Carlo realizations, the results of which are also given 
in Table |4] The test reveals that the correlation of aco with either 
Tdust or /60//100 is stronger than with Alog(SSFR)M5 • 

For the three relations shown in Figure [5]3-d, the constant- 
aco model always leads to higher reduced values than the 
step and linear functions, and can thus be rejected. For the 
(^co-feo/ fioo correlation, the step and linear models are statis- 
tically undistinguishable, while in the case of the aco-Tdust and 
Qfco-Alog(SSFR)M5 correlations, the step function is statistically 
slightly better than the linear model. However, there is a large 
overlap in the reduced x'^ distributions of the two models, pre- 
venting us to adopt or reject either. We note that in the case of 
the step function, the values of a^co below and above the break 
diff'er by a factor of 5, consistent with the relative diff'erence be- 
tween the Galactic and starburst aco values measured in the lo- 
cal Universe (see, e.g., >Solomon et al.|199 7 ). Although we can- 
not distinguish between a step and a linear model, we can unam- 
biguously conclude that aco strongly correlates with Tdust and 
/60//100 and in a weaker way with Alog(SSFR)M5. Before dis- 
cussing the meaning of these correlations (Section [6]), we ffrst 
discuss their robustness against systematic biases in the quanti- 
ties involved in Eq.jT] i.e., ^gdrO^o). ^dust and yUQ. 
(i) Uncertainties on the ^gdr(/^o) relation : 

We have assumed that the z ~ 1 GOODS-N star-forming 
discs follow the local ^gdrO^o) relation derived by Leroy et al^ 
( |2011 ). Instead, this relation might evolve with redshift, both in 
terms of normalization and slope. Because aco Hnearly scales 
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Fig. 5. aco(i-o) a s fun ction of metallicity (i.e., yUo), dust temperature derived from fitting a single modified blackbody function (i.e., 
Tdusu see Section 4^), rest- frame far-infrared color (i.e., /60//100) and distance with respect to the main sequence of star-formation 
(i.e., Alog(SSFR)M5). Filled and opened circles represent our 4 secure and 3 tentative IRAM CO detections, respectively. Our 2 
IRAM CO upper limits are shown by arrows. BzK galaxies and SMGs are shown by filled sq uares and triangles, respectively. The 
shaded area in panel (a) represents the CKco(i-o)-y^o relation and its Icr uncertaint y as found in Genzel et al. ( [2012 ). The shaded area 
in panel (d) represents the Icr scatter around the MS of Rodig hiero et al. ( 2010|). The grey line in panel (c) represents the empirical 
^3^co(i-o)-/6o//ioo relation observed in the local Universe by Gracia-Carpio et al. fSOll ). Red lines in panel (b), (c) and (d) present 
fits to the data points using a constant (dashed-triple-dotted line), a step (dashed-dotted line) and a linear (solid) function. Red lines 
are intentionally limited to the parameter space in which we have direct constraints. 



Table 4. Parameter values of the cKco-^dust, (^co-/6o//ioo and 
(2co-Alog(SSFR)M5 correlations 



Model 




Q^CO-^dust 


Q'CO-/60//l00 


aco-Alog(SSFR)M5 


Spearman's factor 




-0.81 ±0.17 


-0.81 ±0.17 


-0.58 ± 0.26 


Constant'' 


/^red 

a 


5.2 
0.54 + 0.06 


5.2 
0.54 ± 0.06 


6.2 
0.54 ± 0.06 


Step^ 


^red 

a 
b 
c 


3.1 

1.11 ±0.18 
0.37 ± 0.09 
28.7 ± 2.7 


3.4 
0.93 ± 0.24 
0.32 ±0.14 
0.51 ±0.09 


3.4 
1.03 ±0.15 
0.38 ± 0.08 
0.18 ±0.09 


Linear"^ 


^red 

a 
b 


3.5 
2.21 ±0.53 
-0.049 ± 0.017 


3.6 
1.61 ±0.29 
-1.74 ±0.49 


4.2 
0.86 ± 0.09 
-0.57 ±0.14 


Notes. 

^ log(Q'co) = a 
^ log(Q'co) = « if 


x<c 


and log(Q'co) = ^ 


if x>c, where 


•^=7'dust, /60//100 or 



Alog(SSFR)M5 

^ log(aco) ^a + bxx, where x = T^ust, feo/fwo or Alog(SSFR)M5 



with ^gdrO^o) (see Eq. [7]), a diff'erent ^gdrO^o) relation at high- 
redshift would scale upward or downward the inferred values 
of Qfco- Since all but two of our galaxies are situated in a 
narrow range of redshifts (i.e., 1.0 < z < 1.6) and a nar- 
row range of metallicities, a change in the normalization of the 
^gdrO^o) relation with redshift would only translate into a re- 
normalization of our correlations, leaving their slopes mostly 
unchanged. Furthermore, since galaxies with cold and hot dust 
temperatures span similar metallicity ranges (Fig. [5j it would 
also be the case using metallicities inferred with the "funda- 
mental metallicity relation" of Mannucci et al. 2010, see point 
(iii)), a change of the slope of the ^gdr(/^o) should not aff'ect 
significantly our correlations. For example, the correlation of 
Qfco with Tdust, /60//100 and Alog(SSFR)M5 is recovered if us- 
ing a constant value for the gas-to-dust ratio (150) rather than 
the metallicity dependent ^gdrO^o)- To produce the relations of 
Fig. |5j ^gdrO^o) would have to vary diff'erentially for galaxies 
with cold and hot dust temperatures. High dust temperatures are 



usually associated with the dense star-fo rming regions of ma- 
jor merger-driven starburst s ([Sanders etal. 2003 ; Chapman et al.| 
2003 ; Ma gnelli et al.|2012[ [Hay ward et al.|20T2f . In these denser 
environments, the depletion of heavy elements from the gas 
phase is expected to increase ( Jenkins |[2009 ), leading to lower 



^gdrO^o) values. Values of ^gdrO^o) in the dense environment 
of local ULIRGs are inde ed found to be smaller than in nor- 
mal star-forming galaxies ([Solomon et al. 1997[ [Wilson et al. 
[20081 jSantini et al.[[2QT0t [Clements et al. 112010] ); by a factor of 
at most ~ 2. Taking into account this variation of (^gdrO^o) with 
density would however strengthen the correlation of a^co on e.g. 
Tdust- Indeed, Eq.[7]clearly shows that a lower value of ^gdrO^o) 
for galaxies with high densities and hot dust would result in an 
even smaller a^co compared to the normal star- forming galaxies 
with colder dust temperatures. We therefore conclude that it is 
unlikely that the uncertainties on the high-redshift ^gdrO^o) re- 
lation can artificially produce the observed correlation of aco 
with Tdust, /60//100 and Alog(SSFR)M5 • 
( ii) Uncertainties on our dust mass estimates : 

We have followed the approach of Leroy et al. ( 2011[ ) in fol- 
lowing DL07 to infer dust masses. Consequently, any biases in- 
troduced by the DL07 method would aff'ect both our dust mass 
estimates and their ^gdrO^o) prescription. This consistency (also 
in term of rest- frame wavelength coverage, see Section [43| ) pro- 
tects aco from systematic biases introduced by the methodol- 
ogy selected to calculate the dust masses. To erase our correla- 
tions, any dust mass bias would have to be specific to our galaxy 
sample, and correlate w ith their dust tempe ratures. Compared 
to the objects studied in Leroy et al. ( 2011[ ), one specificity of 
our sample is that it also includes galaxies which are analogues 
of local-ULIRG, i.e., merger-driven starbursts like SMGs (see, 
e.g., Tacconi et al . 2008; Magn elli et al. 2012). Because these 
merger-driven starbursts have denser environments, they might 
have high er dust emissivity than normal star-forming galaxies 
(see, e.g., Arce & Goodman 1999VDutr a et al.|2003t Cambresy 
et al. 2005 ; Michalowski et al. 2010 ). Therefore, using the DL07 
method, we may have overestimated the dust masses of our 
merger-driven starbursts with hot dust temperatures. Once again. 
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because aco linearly scales with Mdust, correcting for this pos- 
sible bias will not erase our correlations, but rather strengthen 
them. In any case, we note that Magnelli et al.| ( |2Q1Q| |2Q12| ) 
do not find any significant diff'erences in the dust emissivi ty of 
high-redshift SMGs and lensed-SMGs. |Magdis et aL] ( [2QT2l ) find 
that the absence of millimeter measurements might alfect in a 
diff'erent way the dust mass estimates of normal MS galaxies 
and SMGs. While for normal MS galaxies the dust masses esti- 
mated with or without millimeter measurements are fully con- 
sistent, for SMGs they find that without millimeter measure- 
ments dust masses seem to be overestimated by a factor of ~ 2. 
Among the three SMGs of our sample only one has a millimeter 
measurement. In any case, because this eff'ect would decrease 
the inferred aco value of our SMGs, it would strengthen our 
Q^co-^dust (/60//100) correlation. Finally, one might argue that 
Tdust and Mdust being both estimated from the same FIR mea- 
surements, they are artificially correlated. For a given set of 
FIR measurements, higher Tdust would lead to lower values of 
Mdust and consequently to lower values of aco- However, we re- 
mind the reader that here this eff'ect is at least reduced by the 
fact that Tdust and Mdust are not directly derived from the same 
method, the former being derived using a single modified black- 
body model and the latter being derived using the DL07 model. 
We also note that the aco- feo/f mo and Qfco-AlogCSSFR)^^ cor- 
relations are not aff'ected by this eff'ect. From all these analyses 
we conclude that uncertainties on our dust mass estimates could 
not explain the observed correlations of aco with Tdust. /60//100 
and Alog(SSFR)M5. 

( in) Uncertainties on our metallicity estimates : 

Mannuc ci et al.| ("2010) found that the metallicity of lo- 
cal star-forming galaxies does not only depend on their stellar 
masses but also on their SFRs. Using metallicities inferred with 
the ''fundame ntal metallicity relatio n" of [Mannucci et al.| ( |2010| ) 
and revised in |Genzel et al. ( 2012| ) for z ~ 2 star-forming galax- 
ies, we find no significant change in our correlations. As a sec- 
ond sanity check we assume a constant value for ^gdrO^o) (i-c, 
independent of the metallicity) and verify that correlations are 
still observed. Finally we note that to explain the Ofco-T'dust cor- 
relation with a metallicity eff'ect, all galaxies with cold and hot 
dust temperature should have low and high metallicities, respec- 
tively, and the metallicity range probed by these two galaxy pop- 
ulations should be of ~1 dex. This is unlikely, given the small 
uncertainties on our stellar mass estimates, and the dispersion of 
the mass-metallicity relation. 

We therefore conclude that uncertainties on the ^gdrO^o) re- 
lation, on our measured dust masses or on the inferred metal- 
licities, are unlikely to artificially produce the a^co correlations 
shown in |5] The normalization and slope of the relations are 
however subject to some of these biases. We find a median value 
for aco of 2.5 at the high dust temperature end (rdust > 30 K), 
compared to the value of 1.0 found in local maj or mergers 
dSolomon et"ai[T997l pownes & S olomon] [T9981). Likewise, 
at low dust temperatures (rdust < 30 K) typical of normal star- 
forming galaxies, the median value of aco is 10, while a value 
of 4.34 Mq (K km s~^ pc^)~^ is found in the Mill^ Way ([Strong 



ular phase, we might consider this bias to be spec ific to their 
study. Correcting downward the ^gdrO^o) relation of [Leroy et al. 



& Mattox|1996||Dame et al.|2001t[Leroy et al.|201 ll[Abdo et al. 
20T0I ): 

Other than the sources of uncertainty listed above, we iden- 
tify two other factors that may lead to the systematic overesti- 
mation of aco by a factor of ~2: (i) jLeroy et al.j ( [20TT| ) argue 
that, at fixed metallicity, changes in the dust emissivity and ^qdr 
between the atomic and molecular phases of the ISM may have 
biased upward, by a factor of 1.5 - 2, their estimate of ^gdrO^o)- 
Because the gas in high-redshift galaxies is mainly in a molec- 



(201 1 ) by a factor of 1.5 - 2 would bring our values aco in line 
with those in analogous z = systems, (ii) If indeed the atomic 
gas phase in high-redshift galaxies is negligible, qualitatively its 
inclusion in Eq. [7] would systematically lead to lower values of 
o^co- 



6. Discussion 

Using PdBI C0(2-l) measurements and far-infrared observa- 
tions from HerscheU we estimate the value of aco for a sam- 
ple of 17 high-redshift galaxies, sampling a large region of 
the SFR-M* plane. We find that aco correlates with dust tem- 
peratures (i.e., Tdust) and far-infrared color (i.e., feo/fmo)- A 
weaker correlation with the distance from the main sequence 
of star-formation (i.e., Alog(SSFR)M5) is also observed. None 
of these correlations could be explained by systematic biases 
in our study, although their exact parametrization still remains 
uncertain. Due to these uncertainties, here we only discuss our 
correlations in terms of relative variations of aco with dust tem- 
perature, far-infrared color, and distance of a galaxy from the 
MS. 



6.1. The aco-Tdust stnd aco-feo/fm correlations 

Although low CO rotational lines are almost always very opti- 
cally thick ( Dickman efaT 1986; S olomon et ar[[T987 ), it has 
been shown in the MW and in nearby star-forming galaxies that 
CO line fluxes can be used to derive total molecular gas masses. 
The proportionality factor linking Leo to Mh2, Q^co. can be the- 
oretically inferred using a "virial" approach. At a given metal- 
licity and for a self-gravitating (i.e. virialized) giant molecular 
cloud (GMC), aco scales as: 



aco 



(9) 



where {n(H2)} is the mean hydrogen density and Ti, is the CO 
intrinsic brightness temperature (e.g., [Solomon et al.|1987[ ). 

Assuming that gas and dust are coupled in GMCs (i.e. oc 
^dust), equation [9] may appear to suggest an anti-correlation be- 
tween aco and Tdust- However, observations of clouds in diff'er- 
ent environments of near- solar metallicity (in the MW and in 
other nearby normal star-forming galaxies) suggest that varia- 
tions of {n{H2)) are compensated by temperature changes, such 
that ac o remains constant, with variations of at most a factor of 
2 (e.g . [Strong & Mattox|p^96l [Dame et al.[[2Q0H [Abdo et al 
[20TQI ). 

In the local Universe, the only galaxies with near-solar 
metallicities yet significantly diff'erent values of aco are the 
ULIRGs. In these extreme merger-driven starbursts, the as- 
sumption that molecular gas clouds are in virial equilibrium 
breaks down, and measurements of aco suggest a value reduced 
by a factor of ~ 5 compared to the Galactic value (see, e.g., 
[Solomon et al."19 97l[Downes & Solomon[[1998l[Scoville et al.[ 
1997; Papadop oulos et al.[2012 ). The scaling between ofco, 
(^(7/2)), and T}j however still holds, albeit with a lower propor- 
tionality factor ([Solomon et al.| 19971 [Pownes & Solomon|1998| 



[Sakamoto et al.[[1999| ), indicating that the low aco values are 
mostly caused by a modification of the dynamics of their GMCs 
rather than by the increase in gas/dust temperature (see, e.g.. 



Chapman et a:[[2QQ3l [Clements et ar[[2QTQ| [Wang et al.[[20TT] ). 

These results suggest a step function of aco as a function of 



12 



Magnelli et al.: Molecular gas properties of galaxies in the SFR-M* plane 




Fig. 6. z-band (i.e., F850LP) images from the Hubble Space Telescope for galaxies in our sample with 1.0 < z < 1.6 (i.e., all but two 
galaxies of our sample). Each image is 40kpc x 40kpc, and galaxies are ordered based on their dust temperatures. On average, the 
galaxies with Tdust < 28 K have disk-like morphologies (i.e., PEPJ123709+621507, PEPJ123759+621732, PEPJ123615+621008, 
PEPJ123646+621141, PEPJ123750+621600) while most of the galaxies with Tdust > 33 K are more compact and irregular/clumpy. 



Tdust; normal SFGs dominated by virialized GMCs, which in 
general have cold dust temperatures, have a^co = 4.36 M© (K 
km s"^ pc^)"\ while merger-driven starbursts, which always 
exhibit larger Tdust than normal SFGs, would have ^co ~ 1 
Mq (K km s"^ pc^)"^ This interpretation qualitatively agrees 
with the step model fitted to our ^co-^dust and Q^co-/6o//ioo 
correlations. It also agrees quantitatively with it because in 
this model we find that galaxies with hot dust temperatures 
have aco values five times lower than those with cold dust 
temperatures. Finally, we also note that this interpretation might 
also agree with a simple visual inspection of the morphology 
of our galaxies (Fig. 6]). Galaxies with cold dust temperatures 
(i.e., Tdust < 28 K) have on average a disk-like morphology 
characteristic of normal SFGs (i.e., PEPJ123709+621507, 
PEPJ123759+621732, PEPJ123615-h621008, 
PEPJ123646-F621141, PEPJ123750-h621600), while galax- 
ies with hot dust temperatures tend to be more compact and 
have disturbed morphologies (i.e., Bzk-4lll, Bzk-21000, Bzk- 
17999) consistent with being major-mergers. Of course this 
conclusion has to be taken with caution because of low statistics 
of our sample and the difl&culties encountered by any visual 
morphological classification. 



'2QQ5";"Bo urnaud et al.|2010| ). High-redshift galaxies may thus ex- 
ist through a state where their ISM is a mixture of virialized and 
un-virialized GMCs (this mixture could either be observed at 
all physical scales, or else as a central region with un-virialized 
GMCs surrounded at larger scales by virialized GMCs). The ef- 
fective Tdust and aco derived from our unresolved observations 
would then correspond to a simple linear combination of these 
two "types" of star-forming regions, the cold virialized and hot 
un-virialized GMCs, creating the apparent smooth evolution of 
aco with Tdust (or equivalently feo/fioo)- The existence of two 
"types" of star-forming regions in interacting galaxies has al- 
ready been found locally in the spa tially resolved obser vation s 
of Arp 158 (Boquien et al. [20TTt see also Wei et al. [2QT2|). 
While most star-forming regions of Arp 158 follow closely the 
Schmidt-Kennicutt relation found in spiral galaxies, its nuclear 
starburst and the tip of one of the tidal tails follow more closely 
the Schmidt-Kennicutt relation observed in merger-driven star- 
bursts. The existence of these two "types" of star-forming re- 
gions within the same galaxy seems to favor the interpretation 
where the smooth ^co-^dust relation is due to our spatially unre- 
solved observations. 



Although satisfactory for our step model, this framework 
does not directly explain the possible smooth evolution of a^co 
with Tdust (or equivalently feo/fioo)^ but can be extended to do 
so. High-redshift gal axies hav e larger gas f ractions than in the 
local Unive rse (Tacconi et al.||201Q, ,Daddi eTaLllMTOl |Geach| 



et al.|20lT] ). Under these conditions, simple grav itational insta- 



bilities might locally create un-virialized GMCs ( Springel et al. 



Using hydrodynamic simulations coupled with radiative 
transfer calculations, Narayanan et al. (2012) also find that aco 
might evolve continuously from the Galactic to the merger- 
driven value. They argue that at high redshifts gas rich galaxies 
with their large gravitational instabilities probe a larger dynamic 
range in ISM properties (gas/dust temperatures, density and ve- 
locity dispersion) than local galaxies, leading to a smooth evo- 
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lution of Qfco- In particular, they find that mergers of diff'erent 
mass ratios lead to diff'erent conversion factors. 

While we cannot formally discriminate between a step or 
a smooth evolution of aco with the dust temperature or the 
/60//100 ratio, we can conclude that either quantity can be used 
to select an appropriate value of aco; galaxies with hot dust 
temperatures unambiguously have lower aco than galaxies with 
cold dust temperatures. The main advantage of Tdust for select- 
ing the value of aco is that it does not rely on the morphological 
classification of major-mergers/disky SFGs, which can be diffi- 
cult at high z because clumpy SFGs may be confused with major 
mergers. Moreover, even if the visual classifications of major- 
mergers/disky SFGs would be fully accurate, they would still 
be ineflftcient in distinguishing starbursting-mergers (i.e., with 
likely low aco values) from late mergers with normal mode of 
star-formation (i.e., with likely high Ofco values). Indeed, we ob- 
serve both locally and at high redshift that galaxies with unam- 
biguous morphological signs of merging are not always associ- 
ated with enhanced star formation activity (e.g . |Bushouse et al 



1988t|Dasyra et al.|2QQ6t|Kartaltepe et al.|201 2'). In that case, our 



Tdust criterion would be more indicative of the ISM conditions 
prevailing in these mergers. We discuss further the implications 
of this result in section [631 



6.2. The a CO -^logiSSFK) MS correlation 

Currently, it is commonly interpreted that galaxies situated on 
the main sequence of the SFR-M=, are consistent with a secu- 
lar mode of star formation, while star-forming galaxies located 
far above the MS are powered by strong starbursts with short 
duty-cycles, likely triggered by major mergers (Elbaz et al. 201 1 



Nordon et al.|2012[|Rodighiero et al.|2QTT||Wuyts et al.|2011b 



Magnelli et al.||2012[ Magnelli et al., in prep.). In that view, a 
correlation between aco and Alog(SSFR)M5 can therefore be 
expected: on-MS galaxies should have high aco values while 
above-MS galaxies should have low a^co values as those seen 
in local ULIRGs. Such a correlation is indeed observed in the 
right panel of Fig. |5j though with a lower Spearman's correla- 
tion factor than those found for Tdust or feo/fioo- However, we 
note that this lower Spearman's correlation factor might simply 
reflect the larger uncertainties in deriving Alog(SSFR)M5 than 
^dust or /60//100 (i-G-. uncertainties on M*, SFR and on the exact 
localization of the MS). 

Although we cannot say whether aco depends on 
Alog(SSFR)M5 as a smooth or a step function, we can unam- 
biguously conclude that galaxies situated far above the MS have 
lower value of aco than galaxies situated on the MS. Using 
Alog(SSFR)M5 as a proxy of the aco value will be of great im- 
portance in the near future especially where deep PACS-SPIRE 
Herschel observations are missing to derive r^ust and feo/fioo- 

The observed Qfco-Alog(SSFR)M5 correlation also directly 
confirms the current interpretation that galaxies on the main se- 
quence are in a secular/steadier mode of star-formation, while 
off-MS galaxies are in a s tarbursting mode l ikely t riggered by 
major mergers. However, [Narayanan et al. ( 2Q12| ) argue that 
high-z galaxies, having high gas surface densities and strong 
gravitational instabilities, could have low aco values without re- 
quiring a major-merger. In any case, our observations allow us to 
conclude that the different mechanisms having the power to raise 
a galaxy above the main sequence must also affect the mode of 
star-formation and the physical conditions prevailing in the star- 
forming regions. This interpretation is also supported by the ob- 
servation that galaxies on and above the star formation MS have 
significantly different PAH abundances ( Nordon et al.|2012j ) and 



star-formation rate surface densities ( |Wuyts et al.|2011b"| ), other 
tracers of physical conditions in the ISM. 



6.3. Implications for future work 

Our observations confirm that at high redshift, just like in the 
local Universe, the value of aco varies between galaxies with 
different star formation properties, corro borating results based 
on dynamical mass estimates in SMGs ( [Tacconi et al. 2008). 
Applying a Galactic aco to starbursting galaxies would not only 
lead to unrealistically high gas mass fractions ( Genzel et al.| 
2010 ), but also to high ^gdrO^o) values, inconsistent with local 
observations and chemical evolution models. 

More importantly, our results provide an empiral way to se- 
lect the appropriate value of aco^ using T^ust or feo/fioo- These 
indicators have the great advantage of being less subjective than 
visual morphological classifications of mergers/SFGs, and more 
indicative of ISM conditions than a fixed Ljr criterion. In the 
absence of far-infrared measurements to derive dust properties, 
the off'set of a galaxy from the star formation main sequence, 
Alog(SSFR)M5 , is also a valid indicator. 

The galaxies studied here all have near-solar metallicities, 
as estimated using the mass-metallicity relation. However, both 
observations and models show that a^co increases in low metal- 
licity environments due to the photodissociation of CO, which 
lacks the power of H2 to self-shield against radiation even when 
dust shielding becomes deficient (e.g. Genzel et al. 20 12[ [Glover 



& Mac Low|2QTTt [Narayanan et al.|2011[[Schruba et al.|2012[ ) 

The presence of an orco-rdust correlation has yet to be tested at 
low metallicities, where the main driver of a^co is the availability 
or not of dust shielding. One could expect aco to vary along a 
plane defined by T^ust and metallicity, as suggested by the simu- 
lations of Narayan an et al.| ('2012 ), where a correlation between 
Of CO. gas surface density and merger mass ratio is observed even 
in low metallicity environments. 

The choice of aco also has implications for the Schmidt- 
Kennicutt star formation relation, where SF Gs and starbursting 
major-mergers form two distinct sequences ( [Genzel et al.|2010[ 
Daddi et al.|2010t but see Ivison et al. [2011| ). If aco evolves with 
the dust temperature (or alternatively /eoZ/ioo or Alog(SSFR)M5 ) 
as a step function, it will not modify this observation. Adopting 
a value of aco that varies smoothly with e.g. Tdust would not 
significantly change those results. However, it would lead to a 
smoother transition where starbursting mergers with lower mass 
ratios or gas-rich disks with extreme gravitational instabilities 
would lie between these two sequences. 



7. Summary 

We have combined deep Herschel PACS/SPIRE imaging and 
IRAM PdBI C0(2-l) measurements for a sample of 17 galax- 
ies at z > 1.0, selected out of the GOODS-N field to span a 
broad region of the SFR-M* plane. The far-infrared data were 
use d to d erive accurate dust masses from the model of lDraine"&l 
Li ( 2007| ), while dust temperatures are recovered from the fitting 
of modified blackbody functions of fixed dust emissivity = 1.5. 

2QTT| for the 



Using the prescription of |Leroy et al. 



metallicity-dependence of the gas-to-dust ratio, the measured 
dust masses and the CO luminosities, we inferred the value of 
Of CO in each of the galaxies. The conversion factor is found 
to vary systematically with dust temperature and far-infrared 
color, two tracers of the physical conditions in star-forming re- 
gions. We also found a weaker correlation between aco and 
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Alog(SSFR)M5, the distance from the sequence traced by star- 
forming galaxies in the SFR-M*plane. We have investigated in 
detail possible systematic biases in the analysis, but find none 
that can convincingly explain these correlations. 

The results are consistent with previous observations indicat- 
ing that in normal star-forming galaxies of near-solar metallici- 
ties aco is mostly constant as long as star formation takes place 
in well virialized GMCs, while in major merger-driven starbursts 
where gas is compressed and dust temperatures are higher, the 
value of aco is decreased by a factor of ~ 5. However, our mea- 
surements of aco are consistent not only with a step function, 
but also with being a smooth function of Tdust- This smooth evo- 
lution may come from the fact that we are working with inte- 
grated measurements tracing the average conditions in the galax- 
ies. Either merging events with low mass ratios, or gravitational 
instabilities within the very gas-rich discs could lead to a situa- 
tion where, within a single galaxy, some star formation is taking 
place within normal GMCs while some stars are also formed in 
dense, hot regions. The integrated values of aco and T^ust would 
then be expected to vary anywhere between the Galactic and the 
"starburst" values, producing the smooth relation we observe. 
Alternatively, conditions may be changing at the scale of the star 
forming regions such that the acQ-^'dust relation would be recov- 
ered by resolved observations at the sub-kiloparsec scale. 

Regardless of the exact functional form of the dependence 
of aco on r^ust, our study clearly shows that dust temperature 
(or equivalently the /60//100 ratio) can be used to assign a sin- 
gle conversion factor value to any galaxy. Such a criterion is 
valuable, as it traces the physical conditions of the ISM more 
reliably than a simple SFR or Lir criterion. It is also valuable 
because it is less subjective than methods requiring challenging 
visual classifications of major mergers which furthermore can- 
not easily distinguish between mergers in a starbursting phase 
from those in a normal phase of star-formation. In the absence of 
far-infrared measurements to derive dust properties, the distance 
of a galaxy from the star formation main sequence can be used to 
make a choice between a Galactic conversion factor (aco= 4.36 
Mq (K km s"^ pc^)~\ including the contribution of Helium) and 
a "starburst" value (aco= 1-0 M© (K km s"^ pc^)"^). While most 
galaxies located above the main sequence are major mergers, 
other mechanisms may generate bursts of star formation pow- 
erful enough to displace the galaxies in the SFR-M* plane. Our 
results however indicate that any such mechanism does so by 
modifying the effective mode of star formation and the physical 
conditions prevailing in the sites of star formation. 
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